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Abstract 


The main goal of this study was to investigate the effective factors on the flooding of Gorganrood 
River in Aq’Qala County with emphasis on March 2019 Flood event. In this research 18 sub-criteria 
that play imperious role in creating flood inundation were identified. According to time series flood 
mapping, technical process of satellite data, field survey, and consulting with flood’s experts, the 
sub-criteria were compared in each pairs using the AHP questionnaire, and their final weights were 
determined using Expert Choice software. Topographic and morphological sub-criteria with weights 
of 0.353 and 0.333 had the most impact on the occurrence of inundating, respectively. Based on the 
results, 14.19% of area of Aq’Qala county have located in very low hazard zone (the western and the 
southern regions due to their suitable topography and distance from the Gorganrood river), 18.1% in 
low hazard zone (the northwest of the Aq’Qala due to existence of swampland and salty plains), 
11.95% in moderate hazard zone (the northern of the Aq’Qala due to swampy and hilly nature and 
southern of the Aq’Qala), 29.66% in high hazard zone (the southeast of Aq’ Qala due to closeness to 
the Gorganrood river buffer zones and tributaries of the Qarasou river, and the western areas due to 
the unsuitable morphological conditions and low gradient of Gorganrood bed river), 26.1% in very- 
high hazard zone (the middle of Aq’Qala due to plains’ topography, and qualitative buffer zones of 
the Gorganrood river due to unsuitable morphological conditions). Canalization of outer arch 
meanders in Sallaq Yilqi to Gangirma villages, widening and dredging the river channel at Aq’Qala 
historic bridge, drainage and digging a detour channel according to topographic conditions in middle 
plains can mitigate the flood hazard in the Aq’ Qala County. 
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مقاله پژوهشی 


تعیین عوامل مؤثر بر سیل خیزی رودخانه گرگانرود و ریز پهنه‌بندی خطر سیلاب شهرستان آق‌قلا با 
استفاده از تحلیل سلسله مراتبی AHP‏ 
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حکیده 


هدف از تحقیق حاضر بررسی معیارهای موثر در سیل‌خیزی رودخانه گرگانرود در محدوده شهرستان 
آققلا است. پنج معیار GIS‏ و ۱۸ زیرمعیار فرعی موثر در ایجاد سیلاب رودخانه گرگانرود و شسهرستان 
آق‌قلاه شناسایی و با روش سلسله‌مراتبی AHP‏ معیارها و زیرمعیارها وزن‌دهی شده و لایه‌های آن‌ها در 
نرم‌افزار ArcMap‏ همپوشانی و نقشة خطر سیل شهرستان آق‌قلا تهیه شد. معیارهای توپوگرافی و 
مورفولوژی به ترئیب با اوزان ۰/۳۵۳ و ۰/۲۳۳ بیشترین تأثیر در به وجود آمدن سیلاب را داشتند. ۱۶/۱۹ 
درصد از مساحت شهرستان آققلا در پهنه خطر خیلی کم (مناطق غربی و جنوبی شهرستان با توجه به 
توپوگرافی مناسب و دور بودن از حریم رودخانه)» ۱۸/۱ درصد در digs‏ خطرکم (شمال‌غرب شهرستان با 
توجه به وجود مناطق باتلاقی و شوره‌زار)» ۱۱/۹۵ درصد در پهنه خطرمتوسط (شمال شهرستان با توجه به 
باتلاقی و تپه‌ماهوری بودن منطقه و جنوب شهرستان). ۲۹/۱۱ درصد در پهنه خطرزیاد (جنوب‌شرق 
شهرستان به دلیل Soy‏ به رودخانةٌ گرگانرود و زیرشاخه‌های رودخانه قره‌سو و غرب شهرستان به دلیل 


شرایط نامناسب مورفولوژی و شیب بسیارکم بستر رودخانه) و ۲۷۱ درصد از شهرستان آق‌قلا در پهنه 
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be‏ خیلی‌زیاد (مرکز شهرستان و اطراف رودخانه گرگانرود با توجه به وضعیت توپوگرافی و مورفولوژی 
رودخانه) قرار دارند. دیواره‌سازی و تثبیت قوس خارجی پیچانرودها در محدوده روستاهای سلاق‌یلقی تا 
قانقرمه» لایروبی رودخانه در نزدیکی پل‌هاء تغییر کاربری اراضی متجاوز به حریم رودخانه» لایروبی و 
بسترسازی مناسب در محدوده شهری رودخانه گرگانرود و زهکشی و کانال‌کشی مناسب با شرایط 
توپوگرافی مناطق مرکزی شهرستان از اقداماتی است که به کاهش خطر سیل در شهرستان آق‌قلا کمک 
می‌نماید. 

کلیدواژه‌ها: آق‌قلاه نقشه خطر سیلاب. رودخانه گرگانرود. GIS AHP‏ 


۱-مقدمه 

بدون شک سیلاب به‌عنوان یک بلای طبیعی شناخته شده است که از نظر تلفات جانی و مالی مهیب‌ترین بلای 
طبیعی در جهان محسوب می‌شود (آذرنگ» ۱۳۹۵). کشور ایران با توجه به اقلیم مدیترانه‌ای به لحاظ سیل‌خیزی 
چهارمین کشور دنیاست و سالانه خسارات زیادی را در این زمینه متحمل می‌شود. روند افزایش سیل در پنج دهۀ 
اخیر نشان می‌دهد که تعداد وقوع سیل در دهه ٩۰‏ نسبت به دهه be‏ تقریباً ۱۲ برابر شده است (عبدی‌نژاد. ۱۳۹۶). 
ای ass‏ کی غ ع is ab‏ یشان انم ولا یا مس تام Gad‏ ا همست از 
طریق تخریب پوشش گیاهی» کاربری غیراصولی اراضی: توسعه سطوح غیرقابل نفوذ و دیگر موارد. احتمال 
یل یو راهن اطق کا و را Aes) Saal se‏ ان :۱۳۸۵ زر ی شریان Dis‏ کر dale‏ 
و زمین‌های بیشتری در هنگام طغیان رودخانه‌ها تحت تأثیر قرار می‌گیرند و گاهی حسارت جانی و مالی زیادی را به 
بار می‌آورد (حسین‌زاده, ۱۳۹۶. افزایش سطوح نفوذناپذی که ناشضی از احداث ساختمان در خاک‌های نفوذپذیر 
است. طبعاً از سطوح نفوذپذیر حوضه که قادر به جذب بخشی از بارندگی است. کاسته و درنتیجه بر حجم کل 
رواناب شهر افزوده است و باعث ازدیاد حجم و شدت رواناب و وقوع سیلاب در مناطق پایین‌دست می‌شود 
(صالحیء ۱۳۹۲). سیلاب‌ها نه‌تنها در نواحی برون‌شهری Sh‏ در نواحی شهری نیز آسیب‌های فراوانی به بار می‌آورد. 
شدت این آسیب‌ها گاهی در نواحی شهری بیش از نواحی برون‌شهری است (بدیعی‌زاده ۱۳۹۱). در به وجود آمدن 
سیلاب‌های شهری fle‏ مختلفی وجود دارد که می‌توان به عبور رودخانه‌های حفاظت نشده از Ole‏ شهرهاء تجاوز به 
حریم رودحانه‌ها و از بین بردن مورفولوژی dad‏ تغییرات اقلیمی و بارش با شسدت بالا در مدت‌زمان کم افزایش 
زبری سطوح اشاره کرد (بزرگ‌حداد ۱۳۹۲ نظری‌هاء ۱۳۹۲). نخستین گام در طرح‌های مدیریتی سیلاب داشستن 
نقشه‌های پهنه‌بندی خطر سیل است. این نقشه‌ها آگاهی‌هایی درباره گستردگی و عمق منطقه سیل گیر در دوره 
پاز هان مات رادو اا بیان Car ek gab a‏ خرن OU‏ ا EDE‏ ای 
مالی و جانی سیل را آسان می‌سازد (عسگری» ۱۳۹۶). از مهم‌ترین کاربردهای نقشه‌های پهنه‌بندی حطر سیل می توان 


به تعیین بستر و حریم رودخانه‌ها و تعیین حدود گذرگاه سیل و اراضی سیل‌گیر حاشیه رودخانه. مطالعه و توجیه 


سال دهم تعیین عوامل مؤثر بر سیل‌خیزی رودخانه گرگانرود .... ۷ 


اقتصادی طرح‌های عمرانی با توجه به پهنه سیل در دوره بازگشت‌های متعدد و برآورد خسارت سیل اشاره کرد 
(حیدری» ۱۳۸۵). امروزه با استفاده از سیستم اطلاعات جغرافیایی GIS)‏ و داده‌های ماهواره‌ای سنجش ازدور ازجمله 
مدل رقومی ارتفاعی زمین (DEM)‏ و تصاویر ماهواره‌های راداری و اپتیکی» علاوه بر سرعت محاسبات. می‌توان 
محدوده پخش سیلاب را با دقت بیشتری تعبین نمود (دی بالداسور و ES soll‏ ۲۰۱۲ با توجه به اهمیت موضوع 
OL‏ محققین زیادی به مطالعه این موضوع در داخل و خارج از کشور پرداخته‌اند ازجمله. صالحی و همکاران 
(۱۳۹۲)» پهنه‌بندی خطر سیلاب شهری را با استفاده از نرم‌افزار GIS‏ و فرآیند تحلیل سلسله مراتبی فازی AHP‏ در 
منطقه تهران را با استفاده از سه منبع اصلی cala‏ مدل رقومی ارتفاعی DEM‏ به‌منظور فراهم کردن لایه‌های تجمع 
جریان» شیب و ارتفاع و تصاویر سنجنده TM‏ ماهواره لندست. به‌منظور استخراج اراضی نفوذپذیر و نفوذناپذیر و 
فایل رقومی شبکه آبراهه‌ها و زهکشی آب‌های سطحی در سطح شهر را با روش فرایند سلسله مراتبی فازی (ASL)‏ 
مورد ارزیابی و تحلیل قرار داده که نتایج آن را به‌صورت نقشه‌های پهنه‌بندی خطر سیل در محیط نرم‌افزار 618 ارائه 
کردند. علی و همکاران (۲۰۱۹) با استفاده از روش تصمیم‌گیری چندمعیاره AHP‏ و مدل نرخ بازگشست نقشه‌های 
آسیب‌پذیری و خطر سیلاب را برای ناحیه سونداربان (هندوستان) به دست آوردند. هفت شاحص مهم شامل تراز 
ارتفاعی زمین» درصد شیب. شاخص رطوبت. پراکندگی بارش کاربری زمین و پوشش گیاهی و فاصله از رودخانه 
را در نظر گرفتند. مطالعات آن‌ها نشان داد که فاصله از رودخانه و پراکندگی بارش و کاربری زمین و پوشش زمین از 
مهم‌ترین علل به وجود آورنده سیلاب در این منطقه هستند. چاکرابورتی و موخوپادیای* (۲۰۱۹) نقشه‌های پهنه‌بندی 
خطر سیلاب FRI‏ و نقشه‌های آسیب‌پذیری ]۳۷ شهر کوچبهار (هندوستان) را به کمک روش تصمیم‌گیری 
جندمعیاره AHP‏ و GIS‏ به دست آوردند. شاخص‌های استفاده شده در این تحقیق شامل ترازن شیب سطوح. تراکم 
زهکشی فاصله از رودخانه. حداکثر بارش ویژگی‌های ژئومورفیکی و غیره بودند. بر اساس نتایج آن‌ها ۱ درصد 
۹ درصد و ۳۸/۹۶ درصد شهر به ترتیب در مناطق پرخطر خطر متوسط و کم‌خطر قرار دارند. در دهه‌های اخیر 
استان گلستان شاهد سیل‌های Sop‏ چون سال ۱۳۸۰ در حوضه‌آبریز بالادست رودخانه گرگانرود بوده است که با 
خسارت‌ها ge‏ و مالی سنگینی در شهرستان JUS‏ همراه oy‏ همچنین می‌توان به سیل‌های ۱۳۸۱ و ۱۳۸۶ در 
حوضهآبریز رودخانه گرگانرود اشاره کرد که باعث از بین رفتن ۱۵۰۰ هکتار از اراضی و از بین رفتن ۱۵۰ روستا شد 
(قلی‌نژاده ۱۳۹۱). شهرستان WEST‏ به دلیل عبور رودخانه گرگانرود همواره مورد تهدید سیل‌هایی چون ۸۳۸۰ ۳۸۱ 
۲ ۰۱۳۹۱ ۱۳۹۲ و ۱۳۹۸ واقع شده است. Geb‏ گزارش اداره آب منطقه‌ای استان گلستان ۳۵۰ میلی‌متر بارش فقط 
طی ۶۸ ساعت که معادل 11 درصد بارش سالانه استان گلستان بود در حوضه‌آبریز گرگانرود در فرودین ۱۳۹۸ رخ 


داد. که باعث آبگرفتگی بیش از ۰" درصد شهرستان آق‌قلا شد. لذا شناسایی رفتار رودخانه گرگانرود هنگام وقوع و 
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پس از وقوع سیلاب می‌تواند کمک بزرگی به کنترل سیلاب این رودخانه در محدوده شهرستان آق‌قلا بنماید. هدف 
تحقیق حاضر پهنه‌بندی خطر سیل شهرستان آق‌قلا با توجه به وزن‌دهی و مقایسه زوجی معیارهای مؤثر در به وجود 
آمدن پهنه‌های سیلاب به کمک روش تحلیل سلسله مراتبی AHP‏ است. همچنین نقشه نهایی خطر سیل با پهنه‌های 
سیلاب بهار ۱۳۹۸ که از تصاویر ماهواره‌های سنتینل-۲ و لندست-۸ استخراح شده‌اند و تصاویر پهپاد از زمان واقعه 
سیلاب مقایسه شده است. سپس راهکارهای کاهش خطر سیلاب در مناطق مختلف شهرستان با توجه به وضعیت 
موجود ارائه حواهد شد. نقشه خطر سیلاب شهرستان آق‌قلا جزو مهم‌ترین اسناد هشدار و امدادرسانی در برابر حادثه 
سیلاب است. که می‌توان در اختیار مسئولان ذی‌ربط قرار داد تا از آن به‌عنوان سند چشم‌انداز توسعه شهرستان آق‌قلا 


و Ce pte‏ بحران Oe‏ استفاده نمایند. 


۲- مواد و روش 

۱-۲- منطقه موردمطالعه 

شهر آق‌قلا در استان گلستان و مرکز شهرستان MEST‏ است. که در محدوده مختصات جغرافیایی ۳۷ درجه شمالی 
و ۵6 درجه شرقی قرار داد. شهرستان آق‌قلا در طرفین رودخانه گرگانرود و در شمال شهر گرگان قرار دارد JSS)‏ 
) تقریباً ۵۰ درصد اراضی شهرستان WET‏ به دلیل عبور رودخانه گرگانرود دارای کاربری کشاورزی» کشت گندم و 
جو است. اکثر مناطق شمالی و شمال‌غربی شهرستان آق‌قلا از نوع باتلاقی و شوره‌زار که فاقد پوشش MS‏ و 
درختان مناسب است و همچنین مناطق جنوبی شهرستان آق‌قلا مناطق صنعتی و باغ است که به شهرستان گرگان 
متصل می‌شود. اقلیم شهرستان آق‌قلا مانند منطقه ترکمن‌صحرا بیابانی گرم و مرطوب است. رودخانه گرگانرود از 
مهم‌ترین شاخه‌های سیل‌خیز استان گلستان است که طول OF‏ حدود ۲۵۱ کیلومتر است که از دامنه‌های شمالی 
البرزشرقی و دامنه‌های غربی ارتفاعات استان خراسان‌شمالی سرچشمه می‌گیرد. این رودخانه پس از عبور از شهرهای 
گنبدکاووس و آق‌قلا در غرب خواجه‌نفس با تشکیل دلتای بزرگی به دریای خزر می‌ریزد. حدود AO‏ کیلومتر از 
رودخانه گرگانرود از شسهرستان آق‌قلا و به‌طور Gel‏ حدود ۱۲ کیلومتر آن از محیط شهری آق‌قلا عبور می‌کند 
(گنجی» ۱۳۹۸). متوسط آبدهی سللانه Of‏ حدود ٩۲۰‏ میلیون مترمکعب است (گزارش‌های شرکت آب منطقه‌ای 
استان گلستان» ۱۳۹۵). 
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شکل ۱- محدوده موردمطالعه. شهرستان آق‌قلا و رودخانه گرگانرود 


۲-۲- پهنه‌بندی سیلاب بهار ۱۳۹۸ شهرستان آق‌قلا 

در تحقیق حاضر به‌منظور ارزیابی نقشه خطر سیلاب از تصاویر ماهواره‌ای موجود از سیلاب بهار ۱۳۹۸ استفاده 
گردید. ماهوارۂ ایتیکی ٩6011061-72۸‏ از سری ماهواره‌های Sentinel‏ است که توسط سازمان فضایی اروپا در سال 
۵ به فضا پرتاب شد و در ۱۳ باند در محدوده Hb‏ مرئی» فروسرخ نزدیک و موح‌کوتاه فروسرخ از طیف 
الکترومغناطیسی تصویربرداری می‌کند. قدرت KLE‏ این تصاویر در سه دسته ٠۰‏ ۲۰ 2۰ متر است(صادق J Sg‏ 
۵۹ سه تصویر 56011061-2 از زمان آبگرفتگی در تاریخ‌های ۱۹ و ۲٢‏ فروردین و پنجم اردیبهشت دانلود گردید. 
همچنین یک تصاویر از ماهواره Landsat-8‏ در تاریخ ۱۶ فروردین به‌عنوان مکمل برای اسستخراج پهنه‌های 
سری‌زمانی سیلاب استفاده شد. بعد از اعمال پیش‌پردازش‌های لازم با استفاده از شاخص MNDWI!‏ پهنه‌های 
سری زمانی سیلاب استخراج گردیدند (ژو ٣‏ ۲۰۰۳). نقشه سری زمانی پهنه‌های سیلاب در (شکل ۲) و مساحت 


آن‌ها در (جدول A‏ نشان داده شده است. 


جدول ۱- مساحت پهنه‌های سری‌زمانی سیلاب با استفاده از شاخص by» MNDWI‏ به شهرستان WEST‏ 


مساحت (arrests)‏ تاریخ 
VE 1۱۳‏ فروردین ۱۳۹۸ 
۳۳۱/۳۹ ۷ فروردین ۱۳۹۸ 
\Y¥/08A‏ ۳۹ فروردین ۱۳۹۸ 


1 Sadek and Li 
2 Modified Normalized Water Index, MNDWI: (B3-B11)/(B3+B11) 
3 Xu 
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مساحت (کیلومترمریع) تاریخ 
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شکل ۲- پهنه‌های سری زمانی سیلاب با استفاده از شاخص MNDWI‏ شهرستان آققلا 
۳-۲- ریز پهنه‌بندی خطر سیلاب در محیط GIS‏ 


نقش اساسی داشتند شناسایی گردید. سپس با مطالعات میدانی و مشورت با کارشناسان خبره سیل استان گلستان ۱/۸ 
شاخحص She‏ برگزیده شدند. مراحل شماتیک تهیه ریز پهنه‌بندی حطر سیلاب در (شکل ۳) نشان داده شده و در 
ادامه به‌تفصیل بیان شده است. 


82059 , به زیر معیارها 


هم‌پوشانی لایه‌ها و وزن‌دهی به 
معیارها 


شکل ۳- فلوچارت مراحل agi‏ نقشه‌های خطر سیلاب 
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۱-۳-۲- شناسایی عوامل و معیارهای مؤثر در سیل خیزی رودخانه گرگانرود و شهرستان آق‌قلا 

ابتدا شاخحص‌های موثر در سیل‌خیزی تعریف و شناسایی معیارها به‌صورت دقیق طبقه‌بندی و سپس لایه‌های 
مورد نیاز در بستر نرم‌افزار آماده‌سازی شدند. به‌منظور ایجاد پایگاه داده در نرم‌افزار ArcMap‏ اطلاعات دقیق در مورد 
شاخص‌های مورد نظر از سازمان‌های مرتبط شرکت آب منطقه‌ای, اداره آبخیزداری اداره هواشناسی و اداره کل منابع 
طبیعی استان گلستان جمع‌آوری گردید. زیرمعیارهای شیب جهت شیب تراز یا مدل ارتفاع رقومی, با استفاده از 
نقشه DEM, ALOS‏ با دقت ۱۲/۵ متر استخراج گردیدند؛ زیرمعیارهای بارش دماء Gees‏ با استفاده از داده‌های 
alee E‏ کر GG‏ نی غدل سازی eb‏ کی مها رها ووانات و زک کاک بر اسا 
نقشه طبقه‌بندی خاک منطقه (اداره آبخیزداری استان گلستان) تعیین شد. زیرمعیارهای فاصله از مناطق مسکونی» 
شرایط پل‌هاء فاصله از aos,‏ کاربری اراضیی تجاوز به حریم رودخانه» فاصله از پل‌هاء سینوزیه. عرض رودخانه. 
شرایط کف و دیواره‌های (مورفولوژی رودخانه) و شیب کف رودخانه با بازدیدهای میدانی و نرم‌افزار Google‏ 
0 استخراج و با استفاده از نرمافزار AULOCAD‏ محاسبه شدند (گنجی» ۱۳۹۹). این اطلاعات که شامل لایه‌های 
وکتوری و رستری بودند در نرم‌افزار ArcMap‏ به‌عنوان پایگاه داده ایجاد شد. در ادامه نیاز است تمام لایه‌های موجود 
تبدیل به فرمت رستری شوند تا طبقه‌بندی وزن‌دهی و درنهایت هم‌پوشانی صورت گیرد. زیرمعیارها و معیارها در 
دی ماه هراک )1 lis tO‏ 


تراز 


EE ۳۳5 ew EB eB E) 

(Bre) مت‎ e مه‎ EF (E | 
Y-F. 

1# شمالی شرقی غربی جنوبی مسطح 
۲-۶ 
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۲-۳-۲- تهیه پرسشنامه AHP‏ و تعیین وزن معیارها و زیرمعیارها 

یکی از روش‌های ارزیابی و وزن‌دهی به معیارها و زیرمعیارها؛ روش تحلیل سلسهه مراتبی AHP‏ است که 
به‌وسیله ساعتی در سال ۱۹۸۰ مطرح شد. نرم‌افزار ExpertChoice‏ این امکان را فراهم می‌سازد تا تصمیمات مستدل 
و منطقی جهت انتخاب مطلوب برای تصمیمگیرندگان فراهم شود (آقایی» (VAT‏ به‌منظور بررسی اهمیت و درجه 
تأثیر هر یک از معیارها و زیرمعیارها در به وجود آمدن سیلاب در شهرستان آق‌قلا از فرایند تحلیل سلسله مراتبی 
AHP‏ استفاده گردید. بدین منظور با استفاده از طراحی پرسشنامه AHP‏ برای فاکتورهای مورد نظر در مقیاس ۱ تا ٩‏ 
بر حسب میزان اهمیت و مقایسه دودویی و ضریب اهمیت هر یک از فاکتورهای مورد بررسی در اختیار ۱۰ 
eee‏ و کار EL‏ خی ebony Dbl Ol CN‏ مربوطه فراز وا مه ما کین بط را 
کارشناسان به‌منظور وزن‌دهی در نرم‌افزار ExpertChoice‏ وارد گردید. وزن‌های دست آمده بر اساس مقایسه 
دودویی و زوجی هر یک از معیارها با یکدیگر بر اساس نظر کارشناسان است. پس از به دست آوردن وزن نهایی هر 
معیار» بايد به محاسبه نرخ ناسازگاری (ConsistencyRatio)‏ پرداخت. نرخ ناسا زگاری به معنی میزان ارتباط منطقی 
بین ایده‌های مختلف با یکدیگر با توجه به یک ملاک معین است. درصورتی که نرخ ناسازگاری کوچک‌تر از ۰/۱ 
باشد (CR<0.1)‏ سازگاری وجود دار و در غیر این صورت Lb‏ در مقدار ارجحیت‌ها تجدید نظر شود (اصغرپور 
۷ مقدار نرخ ناسازگاری از (رابطه ۱) و (رابطه ۲) به دست می‌آید. 


CI 
— AMax Mn 
CI = a )۲ (رابطه‎ 


IRI (میانگین عناصر بردار سازگاری)-‎ Amax (شاخص ناسا زگاری)- 0 (بردار نرمال‌سازی)-‎ CI ناسا زگاری)-‎ cy) CR 


(شاخص ناسا زگاری تصادفی (جدول ۲)). 


حدول ۲- حدول شاخص ناسازگاری تصادفی 
n ۱ ۲ ۳ é ۵ 1 ۷ ۸ 4 ۱۰‏ 


IRI ۰ ۰ ۰/0۸ ۰/۹۰ \/\Y 24 VYY WEN \/£0 \/0\ 


در صوری که بعد از محاسبات فوق مقدار نرخ ناسازگاری ۰/۱ CRS‏ به دست آمد. نشان دهنده این است که 
ارجحیت‌های وارد شده صحیح است (عطایی OTM‏ بعدازاینکه مقایسات زوجی توسط کارشناسان صورت گرفت؛ 
با جمع‌بندی و میانگین گیری از ارجحیت‌های انجام شده پرسشنامه نهایی براساس مقایسات نهایی وارد نرم‌افزار 
6 گردید. اوزان نهایی به دست آمده معیارها و زیرمعیارها و گزینه‌ها نسبت به همدیگر و نسبت به 
سطح بالاتر خود و درنهایت نسبت به هدف نهایی به دست آمد. اوزان به دست آمده و همچنین نرخ ناسازگاری 


هرکدام در (جدول ۳( نشان داده شده است. 
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جدول ۳- اوزان نهایی به دست آمده با استفاده از نرم‌افزار ExpertChoice‏ 


نرخ ناسازگاری وزن نهایی معیار؟:زیرمعیار وزن معیار و زیرمعیار معیارها و زیرمعیارها 
“or : >‏ توپوگرافی 
vy‏ ۰/۹ شیب 
۰/0/۱ ۰۱۹ ۰/۳۱ جهت شیب 
Vt‏ ۹۳ تراز- ارتفاع 
ا e‏ 2 هیدرواقلیم 
۱ 4 بارش 
LAA ۹‏ دما 
۱/< 
۹ ۹۸ تبخیر 
eto‏ ۳۹ رواناب 
z <‏ عوامل انسانی 
)+/+ 1/۹ فاصله از مناطق مسکونی 
UYN 0‏ شرایط پل‌ها 
YA ۵‏ فاصله از رودخانه 
۱ 
aS‏ ۹۹ کاربری اراضی 
۱ ۰۱۳۵ تجاوز به حریم رودخانه 
“AY ۲‏ فاصله از پل‌ها 
۱/۰۱ = ۱/۹ زمین‌شناسی 
)/< 1۹ ۱ زهکشی خاک 
‘NYT ۰/۰/۱‏ مورفولوژی 
ANI‏ ۰/۳۵۱ سینوزیته رودخانه 
۳ 71۲ عرض رودخانه 
۱ 
A 4‏ شرایط کف و دیواره رودخانه 
۸۲ ۳:۷ شیب بستر رودخانه 


۳-۳-۲- ایجاد پایگاه داده در ArcMap‏ و تبدیل داده‌ها به فرمت رستری 


داده‌های وارده شده به نرم‌افزار ۵۲٩13۲‏ از نوع وکتوری و رستری بوده که برای طبقه‌بندی تصاویر نیاز است 
تمام داده‌ها به فرمت رستری تبدیل شوند. در فرمت رستری هرکدام از گزینه‌های موجود در لایه‌ها به پیکسل تبدیل 
شده که هر کدم از پیکسل‌ها دارای ارزش است. برای طبقه‌بندی نیاز است تا بر اساس ارزش هرکدام از گزینه‌ها یک 
ارزش به آن‌ها داده شود تا درنهایت تصاویر خحروجی» بر اساس ارزش پیکسل سیمبولوژی شوند که نشان دهنده 
اهمیت گزینه مورد نظر است. از میان hay‏ استفاده شده در پژوهش. shay‏ شیب. جهت شیب و تراز جون 


مستقیم از فایل DEM‏ استخراج می گردند دارای فرمت رستری می‌باشند؛ اما داده‌های بارش دماء تبخیر» فاصله از 
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مناطق مسکونی. فاصله از پل‌هاء تجاوز به حریم رودخانه» سینوزیته» عرض SUS‏ شرایط کف و دیواره و شیب کف 
رودخانه به دلیل اینکه از نوع نقطه‌ای هستند و تمام اطلاعات آن‌ها > Attribute Table‏ ذخیره شدهاند با استفاده از 
دستور IDW!‏ برای تبدیل آن‌ها به فرمت رستری استفاده گردید (پرادهان2009). برای تبدیل SY‏ خطی فاصله از 
رودخانه» از دستور فاصله اقلیدسی" و لایه‌های پلی‌گون رواناب کاربری اراضی و زهکشی خاک به فرمت رستری از 
دستور طبقه‌بندی مجدد" استفاده گردید. 

۳-۲-ع- هم پوشانی لایه‌ها و وزن‌دهی به معیارها 

بعدازاینکه تمام لایه‌ها به ده کلاس تقسیم‌بندی گردید. طبقه‌بندی مجدد لایه‌ها به پایان می‌رسد. در طبقه‌بندی» 
کلاس‌یک کمترین ارزش در خطر سیلاب و کلاس‌ده بیشترین ارزش در خطر سیلاب در شهرستان WEST‏ را دارد. 
سپس با استفاده از دستور Weighted Sum‏ لایه‌ها هم‌پو شانی شده و وزن‌های نهایی (جدول ۳) به دست آمده 


مربوط به هر معیار و زیرمعیار به لایه مربوطه داده می‌شود. 


۳- نتایج و cou‏ 

۱1-۳- نقشه ole‏ زیرمعیارهای طبقه‌بندی شده در محیط GIS‏ 

بعد از تبدیل همه لایه‌ها به فرمت رستری نیاز است تا متغیرهای همه لایه‌ها یکسان‌سازی و دارای ارزش شوند. 
ides‏ منظور تمام لایه‌هاء پس از وزن‌دهی و ارزش‌دهی به زیرلایه‌هاء به ده کلاس تقسیم‌بندی گردیدند. کلاس یک 


کمترین ارزش در به وجود آمدن سیلاب و کلاس ده بیشترین ارزش در به وجود آمدن سیلاب در شهرستان آققلا 
را دارند. ada‏ طبقه‌بندی شده زیرمعیارها به تفکیک در (شکل ۵ آورده شده است. 


1 Inverse Distance Weighted 
2 Pradhan 

3 Euclidean Distance 

4 Reclassify 
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شکل 0- )\-\( ay‏ درصد شیب (۲-۱) Ce ay‏ شیب. (۲-۱) ay‏ مدل رقومی ارتفاع؛ (۱-۲ لا یه بارش. 
a ۲-۲(‏ دما )9 a G1) sur a‏ روافاب:۱-۳(۰ لاه فاصله از مناطن مسکوتی: OY)‏ لاه شراط 


پل‌هاء (۳-۳) a‏ فاصله از رودخانه. (4-۳) لایه کاربری اراضی. )0-1( لایه فاصله از پل‌ها؛ (۱-۶) لایه بافت 


خاک؛ (۱-۵) لایه ضریب خمید گی (سینوزیته)؛ (۲-۵) aY‏ عرض a )۳-۵( SUIS‏ شرایط کف و دیواره. (۵- 


)٤‏ لایه درصد شیب بستر رودخانه 
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۲-۳- نقشة نهایی پهنه‌بندی خطر سیلاب در محیط GIS‏ 

در گام آخر بعد از وزن‌دهی معیارها و زیرمعیارهاء لایه‌های تهیه شده به ده کلاس طبقه‌بندی شدند که ارزش ده 
بیشترین خطر سیلاب و ارزش یک کمترین خطر سیلاب را داشت. لایه‌ها با استفاده از ابزار مجموع وزنی 
(Weighted Sum)‏ هم‌پوشانی شدند و وزن هر زیرمعیار به لایه آن اعمال گردید. (شکل A‏ نقشه نهایی خطر 


سیلاب شهرستان و شهر آق‌قلا را نشان می‌دهد. 
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شکل *- EA‏ نهایی خطر سیلاب شهرستان و شهر آق‌قلا 


با توجه به (شسکل )٩‏ مناطق قرمزرنگ و آبی‌رنگ به ترتیب مناطق پرخطر و کم خطر در برابر مسیلاب می‌باشند. 
اکثر مناطق اطراف رودخانه گرگانرود با توجه به معیارهای مورفولوژی و عوامل انسانی در Hg‏ پرخطر قرار گرفتند. 
به دلیل اینکه الگوی رودخانه گرگانرود در محدوده روستای سلاقیلقی تا روستای قانقرمه از نوع پیچانرودی شدید 
است. در سیلاب بهار ۹۸ با توجه به عدم لایروبی رودخانه و فرس‌ایش کناره‌های رودخانه بستر و عرض مفید 
رودخانه کاهش adh‏ و سیلاب جاری در رودخانه در این ناحیه از مدخل رودخانه خارج شده و با توجه به شیب 
شرقی- غربی منطقه از دشت‌های شمالی شهر آق‌قلا به محبط شهری شهر آق‌قلا وارد گردیده است. لذا این منطقه از 
شهرستان در محدوده خطر خبلی‌زیاد قرار می گیرد. همچنین دشت‌های میانی شهرستان با توجه به معیار توپوگرافی و 
eg‏ اراک در مدرد ee‏ قزار ک رتوا ماد رشان دو شبات هار ۱۳۹۸ a‏ ول Side Ke‏ 
آب‌بندان پلمه‌خندان آب موجود در آن به دشت‌های مرکزی شهرستان جریان یافته و با توجه به شیب بسیار پایین 
منطقه (دو در ده‌هزار) سیلاب به‌صورت راکد مانده و منطقه‌ای پرخطر در برابر سیلاب به وجود می‌آید. مناطق شمالی 
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و شمالی‌شرقی شهرستان با توجه به زیرمعیارهای شیب و جهت شیب. باعث عدم > CS‏ سیلاب و وجود مناطق 
باتلاقی باعث به دام افتادن سیلاب خارج شده از رودخانه اترک در این ناحیه می‌شسود لذا باوجوداینکه این منطقه از 
شهرستان فاقد مناطق مسکونی و زمین‌های زراعی است. در پهنه با خطر متوسط قرار گرفت. مناطق غربی و شمال- 
غربی شهرستان با توجه به وجود دشت‌های شوره‌زار و باتلاقی و دور بودن از رودخانه گرگانرود و نبود منازل 
مسکونی در محدوده خطر کم تا خیلی‌کم قرار گرفت. این ناحیه در سیلاب بهار ۱۳۹۸ به Jo‏ طغیان رودخانه اترک 
دچار پهنه‌های سیلابی شده است. مناطق غربی و جنوبی شهرستان به دلیل دور بودن از رودخانه گرگانرود و شیب 
نسبتاً مناسب در محدوده خطر خیلی‌کم قرار گرفتند. 

۳-۳- ارزیایی صحت نقشه خطر سیلاب 

یکی از گام‌های اساسی در ارزیابی صحت نقشه‌های تولید شده توسط سیستم اطلاعات جغرافیایی» مقایسه آن‌ها 
با تصاویر ماهواره‌ای در دسترس و تصاویر پهپادی موجود از زمان واقعه سیلاب است. که جزو دقیق‌ترین مراجع 
جهت صحت نتایج می‌باشند. بدین منظور نقشه خطر سیلاب (شکل 1) با پهنه‌های سری‌زمانی سیلاب استخراج شده 
از تصاویر موجود byl pale‏ ستتینل-۲ و لندست-۸ (شکل (A‏ و تصاویر پهپاد زمان سیلاب (شکل 4( به‌صورت کیفی 
مقایسه گردید. 

بعد از استخراج پهنه‌های سیلاب به کمک تصاویر ماهواره‌ای و همچنین تعیین وزن معیارهای تحقیق بیشترین 
اوزان مربوط به معیارهای توپوگرافی و مورفولوژی به ترتیب ۰/۳۵۳ و ۰/۳۳۳ بودند. با بررسی‌های انجام شسله بر 
روی نقشه پهنه‌بندی خطر سیلاب (شکل )١‏ و پهنه‌های سیلاب استخراج شده از تصاویر ماهواره‌ای (شکل AY‏ 
مشخص شد در مناطقی که معیارهای توپوگرافی و مورفولوژی رودخانه در وضعیت حساسی قرار دارند. باعث ایجاد 
پهنه‌های سیلاب شده است. این مناطق به‌طور حاص در شعاع پنج کیلومتری شهر آققلا در محدوده روستای 
سلاقیلقی تا روستای دوگونچی قرار دارند. که باعث شده است این محدوده در jee‏ خیلی‌زیاد سیلاب قرار گیرد. 
در ادامه به چند عامل مهم در ایجاد پهنه‌های سیلاب در این محدوده اشاره شده است. 

ضریب سینوزیته رودخانه در محدوده روستای سلاق‌یلقی (بازه اول. شکل ۷ برابر ۳/۸۸ است که رودخانه از 
نوع پیچانرودی شدید اسست. اولین پیچانرود در ijk‏ اول با شعاع حدوداً ۸۰ متر باعث خروج سیلاب از مدخل 
رودخانه شده که با توجه به شیب و توپوگرافی و نبود پوشش گیاهی جاذب آب. سیلاب به دشت‌های واقع در شمال 
شهر آق‌قلا جریان پیدا کرده و از مناطق شمالی شهر وارد محیط مرکزی شهر شده است. ضریب سینوزیته رودخانه 
در محدوده روستای قانقرمه (بازه دوم-شسکل ‏ برابر ۲/۰۳۹ است که رودخانه از نوع پیچانرودی شدید است؛ اما با 
توجه به اینکه حجم bles‏ سیلاب در ail‏ اول از مدخل رودخانه خارج شده حجم سیلاب کمتری به EJL‏ دوم 


Ob >‏ پیدا کرده است. در روستای محمدآلق» تصرف به حریم کمی رودخانه باعث تخریب دیواره‌های رودخانه در 
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بعضی از مناطق شده و وجود پل عبور و مرور با تراز سطح عرشه پایین و عدم وجود پایه‌های مناسب مانند مانعی بر 
سر راه جریان باعث بالا آمدن سطح جریان و خروج سیلاب از مدخل شده است. ضریب سینوزیته رودخانه در 
محدوده شهری آق‌قلا (بازه سوم شکل ۷ برابر ۱/۱۶ می‌باشد که رودخانه از نوع پیچانرودی است. اکثر مسیر 
رودخانه در این بازه مستقیم بوده و میانگین شعاع پیچانرودها هشت درصد از ميانگین شعاع پیچانرودهای بازة اول 
بیشتر است. میانگین عرض بستر کانال رودخانه در بازه سوم ۲1/ بیشتر از میانگین عرض بستر کانال رودخانه در 
بازه‌های اول و دوم است. ضریب سینوزیته رودخانه در محدوده روستای دوگونچی (بازه چهارم» شکل V‏ برابر 
۶ که رودخانه از نوع پیچانرودی اسست. میانگین عرض پیچانرودها در این بازه ۳۰/ نسبت به بازه سوم PS‏ 
است. که نشان می‌دهد عرض رودخانه در پیچانرودها کمتر از سایر بازه‌ها است. در این بازه پهنه‌های سیلاب فقط در 
نزدیکی پیچانرودها از مدخل رودخانه خارج شده است. تغیبرات SUSE‏ مقطع رودخانه و تغییر شیب از (سه در ده 
هزار به دو در ده هزار) باعث جدا شدن خطوط جریان شده و تخریب مخصوصاً در کناره‌های پل حوالی روستای 


آق‌تکه‌خان به طرف پایین‌دست رودخانه گرگانرود می‌شود. 


~ 
مناطق مسکونی 
مرکز شهر WEE!‏ 
" رودخانه گرگانرود 


شهر آق‌قلا 


شکل ۷- در حدفاصل cling,‏ سلاق یلقی تا روستای دوگونجی. محدوده خطر خیلی‌زیاد 
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شکل ۸- مقایسه نقشه خطر سیلاب با پهنه‌های سیلاب استخراج شده به کمک تصاویر ماهواره سنتیئل -۲ و 
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شکل -٩‏ مقایسه پهنه‌های پرخطر با تصاویر پهیاد. سیلاب بهار ۱۳۹۸ شهرستان Wal‏ 


£- نتیجه‌گیری و پيشنهادها 
۱-6- پیشنهادهای کاهش آسیب‌پذیری خطر سیلاب 


به‌منظور تفکیک شهرستان آق‌قلا با توجه به خطر سیلاب نقشه (شکل ۱۰) به پنج بازه» حطر خیلی کې حطر کې 
خطر متوسط. خطر زیاد و خطر خیلی‌زیاد تقسیم‌بندی گردید. با توجه به (شکل ۰ مناطق اطراف رودخانه 
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گرگانرود در محدوده سد وشمگیر تا شهر آق‌قلا تا شعاع ده کیلومتر در پهنه خطر خیلی‌زیاد و مناطق شمالی و جنوبی 
شهرستان در پهنه خطر obj‏ قرار دارند. که جزو مناطق حساس در برابر سیلاب محسوب می‌شوند. در ادامه با توجه 
le oh‏ تس هط میات ره رها کون سای واه شاه olen eel‏ درل کی 
پهنه‌های خطر در (جدول )٤‏ نشان داده شده است. 

پهنه با حطر خیلی زیاد (رنگ قرمز): حدود ۲ درصد شهرستان آق‌قلا در محدوده خطر سیلاب قرار گرفت. با 
یشان توش کر یاه SUS‏ کرت مر ی e‏ 
به سمت دشت‌های شمالی که شوره‌زار و باتلاقی است. می‌تواند هنگام سیل تا حد زیادی دبی پیک خروجی از سد 
وشمگیر را به آن مناطق منحرف کند تا حجم سیلاب کمتری وارد رودخانه گرگانرود شود با توجه به ماهیت 
ays‏ کرافت متطفهه زاشکاز ا ا وم AG bese AKA HOG‏ ماند. CSCI‏ 
و هشدار به‌موقع» هنگام سیل با آن به تعامل پردانعت. مناطق مرکزی شهرستان به دلیل وجود آب‌بندان‌های زیاد 
ازجمله آب‌بندان یلمه‌عندان هنگام واقعه سیلاب دبی زیادی از سر ریز غربی سد وشمگیر وارد این آب‌بندان‌ها 
می‌شود. کنترل حجم دبی ورودی به این آب‌بندان‌هاء تثبیت بدنه آن‌ها و زهکشی و حفر کانال‌های جمع‌آوری رواناب 
با توجه به شرایط توپوگرافی موجود راه‌حل مناسبی است. همچنین جابه‌جایی مسیر ریل راه‌آهن یا ایجاد کالورت‌های 
عبور سیلاب در زیر حطوط ریل راه‌آهن باعث جریان سیلاب به دشت‌های شمالی شده که فاقد مناطق مسکونی و 
زمین‌های زراعی است. 

پهنه با حطر زیاد (رنگ زرد): حدود ۳۰ درصد شهرستان آق‌قلا در محدوده خطر سیلاب قرار گرفت. این پهنه‌ها 
در مناطق شمال-شرقی. جنوب-شرقی و مناطق غربی شهرستان WET‏ قرار دارند. در مناطق جنوبی سد وشمگیر به 
علت تجاوز شهرها و روستاها به حریم رودخانه» وجود منازل مسکونی موجود در حریم کمی (۲۰ متر) رودخانه 
گرگانرود هنگام سیل با دوره بازگشت ۵۰ سال به بالا دچار تخریب می‌شوند. جابه‌جایی این منازل و تغییر کاربری 
آن‌ها به زمین کشاورزی و باغ می‌تواند باعث افزایش نفوذپذیری و افزایش ضریب زبری شود همچنین اطلاع‌رسانی 
به‌موقع» سیستم هشدار سیل سد وشمگیر می‌تواند کمک بزرگی به کاهش خطر جانی و مالی کند. مناطق شمال- 
شرقی شهرستان اکثر مناطق شوره‌زار و باتلاقی است که هنگام طغیان رودخانه اترک دچار آبگرفتگی می‌شوند و با 
توجه وجود تیه‌ماهوری منطقه و عدم زهکشی منطقه جریان سیلاب به سمت مناطق غربی و دریای خزر وجود ندارد 
لذا با توجه به اهمیت کمتر این منطقه با توجه به ساير مناطق شهرستان راهکارهای کاهش خطر سیلاب در منطقه 
صرفه اقتصادی نخواهد داشت. 

پهنه با be‏ متوسط (رنگ سبز): حدود ۱۲ درصد شهرستان آق‌قلا در محدوده خطر سیلاب قرار گرفت. مناطق 
جنوبی شهر آق‌قلا به دلیل عبور زیرشاخه‌های رودخانه قره‌سو لایروبی این رودخانه و بهسازی پل‌های بتنی می‌تواند 
کمک بزرگی به کنترل سیل در این نواحی AS‏ نواحی غربی سد وشمگیر به دلیل کانالیزه کردن رودخانه و شیب 
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مناسب در سیلاب با دوره بازگشت ۵۰ ساله دچار آبگرفتگی نشدند. دشت‌های شمالی شهرستان با توجه به اینکه در 
حریم رودخانه اترک قرار دارد. در سیلاب بهار ۹۸ با بیرون‌زدگی سیلاب از SUIS‏ رودخانه اترک» رواناب آن با توجه 
به ماهیت توپوگرفی tlie‏ وارد ناطق شمالی شسهرستان می‌شسود؛ این مناطق به دلیل وجود دشست‌های باتلاقی و 
شوره‌زار مکان مناسبی برای هدایت سیلاب می‌باشند. 

پهنه با خطر کم (رنک آبی): حدود ۱۸ درصد شهرستان آققلا در محدوده خطر سیلاب قرار گرفت. مناطق 
شمال-غربی شهرستان اغلب شامل شوره‌زار و مناطق باتلاقی می‌باشند که کاربری خاصی ندارند و با ایجاد JUS‏ 
انحرافی و زهکشی مناسب حوضه‌آبریز حجم عمده‌ای از رواناب‌های جاری را به سمت دریای خزر هدایت کرد. 

پهنه با حطر خیلی کم (رنگ فیروزه‌ای): حدود ۱۶ درصد شهرستان آق‌قلا در محدوده خطر سیلاب قرار گرفت. 
این مناطق به دلیل اینکه در فاصله دورتری از رودخانه گرگانرود قرار دارند و اغلب کاربری اراضی آن باغ و زمین‌های 
زراعی است و ممکن است در سیل‌های با دوره بازگشت ۵۰ ساله به بالا دچار آبگرفتی خفیف شوند. لذا با انجام 
cle ol‏ ارائه شده این مناطق جزو مناطق ایمن در برابر سیلاب می‌باشند. 


حدول ۶- مساحت و درصد پهنه‌های خطر سیلاب. شهرستان آققلا 


پهنه‌های خطر مساحت (کیلومترمربع) درصد 
پهنه با خطر Yvi A0۵0 Lj le‏ 
پهنه با خطر زياد 299۳/۱۸ ۳۹/۳۹۹ 
پهنه با خطر متوسط ۳۳۳۸/۳/5۸ ۱۱/۳۹۵ 
پهنه با خطر کم ۱۳۳/۹۶ ۱۸/۱ 
پهنه با خطر خیلی کم YUL/0VE‏ ۱/۹ 
مجموع ۸7714 ۱۰۰ 
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شکل ۱۰- نقشه تفکیک شده خطر سیلاب شهرستان آققلا 
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0 جمع‌بندی 

بر اساس uli‏ رودخانه گرگانرود در حدفاصل روستای سلاق‌بلقی تا روستای قانقرمه جزو مناطق حساس از 
نظر مورفولوژی است. در این بازه شاخص سینوزیته برابر ۳/۸۸ تعداد پیچانرود و میانگین شعاع پیچانرودها به ترتیب 
برابر ۱۸ و ٩۸متر‏ اسست (گنجی» ۱۳۹۹) که نسبت به سایر بازه‌ها تعداد پیچانرودها و ميانگین شعاع پیچانرودها به 
ترتیب ۳۰ درصد بیشتر و پنج درصد کمتر است. با بازدیدهای میدانی و مشورت با کارشناسان خبره سیلاب. پنج 
معیار کلی و ۱۸ زیرمعیار موثر در ایجاد سیلاب شهرستان آق‌قلا و رودخانه گرگانرود شسناسایی گردید. بعد از 
وزن‌دهی آن‌ها با روش سلسله مراتبی AHP‏ معیارهای توپوگرافی و عوامل مورفولوژی به ترتیب ۰/۳۵۳ و ۰/۳۳۳ 
بیشترین اوزان را به خود احتصاص دادند. براساس نقشه آسیب‌پذیری خطر سیلاب. مشخص شده حدفاصل سد 
وشمگیر تا ابتدای شهر WET‏ و دشست‌های مرکزی شهرستان به علت شرایط نامناسب توپوگرافی و مورفولوژی 
رودخانه گرگانرود جزو پهنه‌های پرخطر در برابر سیل می‌باششند. با توجه به تقشه تفکیک شده آسیب‌پذیری خطر 
سیلاب مشخص شده ۱2/۱۹ درصد از مساحت شهرستان آق‌قلا در پهنه با خطر خیلی کم و ۱۸/۱ درصد در پهنه با 
خطر کم و ۱۱/۹۵ درصد در پهنه با خطر متوسط و ۲۹/۱ درصد در پهنه با خطر زیاد و درنهایت ۲۷/۱ درصد 


شهرستان در پهنه با خطر خیلیزیاد قرار دارد. 


-A‏ تشکر و قدردانی 
نویسندگان مقاله از اداره منابع طبیعی و آبخیزداری و شرکت آب منطقه‌ای استان گلستان Cob‏ همکاری در تکمیل 
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